Scaling Solutions of Inelastic Boltzmann Equations 
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C/j ' This paper deals with solutions of the nonlinear Boltzmann equation for spatially uniform freely 

cooling inelastic Maxwell models for large times and for large velocities, and the nonuniform con- 
vergence to these limits. We demonstrate how the velocity distribution approaches in the scaling 
limit to a similarity solution with a power law tail for general classes of initial conditions and derive 
a transcendental equation from which the exponents in the tails can be calculated. Moreover on the 
basis of the available analytic and numerical results for inelastic hard spheres and inelastic Maxwell 
models we formulate a conjecture on the approach of the velocity distribution function to a scaling 
form. 



CO . I. INTRODUCTION 

> 

In recent times overpopulation of high energy times in velocity distributions of freely cooling or in driven granular 
fluids has become a focus of attention in laboratory experiments, kinetic theory and computer simulations Jy. In 
kinetic theory, granular fluids out of equilibrium were originally modeled by inelastic hard spheres (IHS) (^J^l, which 
is the proto-typical model for dissipative short ranged hard core interactions. In general, similarity solutions are 
of interest, because they play an important role as asymptotic or limiting solutions of the Boltzmann equation at 
large times or at large velocities, and they frequently show overpopulated high energy tails when compared to the 
omni-present Gaussians. 

Overpopulated tails were first found theoretically by studying scaling or similarity solutions of the nonlinear Enskog- 
Boltzmann equation for the IHS fluid, both in freely evolving IHS systems without energy input Q, as well as in 
driven or fluidized systems and confirmed afterwards by Monte Carlo simulations of the Boltzmann equation 

jj],U j and by laboratory experiments . Overpopulated tails in free IHS fluids || and driven ones have also been 
studied by molecular dynamics simulations of inelastic hard spheres. The observed overpopulations in IHS systems 
are mainly stretched exponentials exp[— Ac b ] with 6=1 m systems without energy input, and b — 3/2 

in driven systems, but for some forms of driving ^,|| 6 = 2 has been observed. For hard sphere systems there is in 
general good agreement between the analytic predictions and numerical or Monte Carlo solutions of the nonlinear 
Boltzmann equation. 



5— I 

About two years ago Boltzmann equations for inelastic Maxwell models have been introduced for the free case 



without energy input in Ref. |12yi3j ] , and for the driven case in Ref. 
received a great deal of attention, because general classes of such so 
unphysical. 



. In these studies similarity solutions have 
utions turned out to be non-positive, and hence 
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The interest in overpopulated tails in elastic Maxwell models has already a long history, and originated from the 
discovery of an exact positive similarity solution of the nonlinear Boltzmann equation for Maxwell molecules, the 



so-call BKW mode 1 17- 20 1 , named after Bobylev, Krook and Wu. The most recent interest in similarity solutions 
of the Boltzmann equation for inelastic Maxwell models (IMM) was also stimulated by the discovery of an exact 
similarity solution |21| for a freely cooling one-dimensional IMM of the form /(c) = (2/w)l/[l + c 2 ] 2 . It is positive, 



has finite mean energy (c 2 ), and an algebraic high energy tail 1/c 4 . The same authors also obtained |2l],|2^] Monte 
Carlo solutions f(y, t) of the nonlinear Boltzmann equation for freely evolving IMM systems with general initial 
distributions both in one and two dimensions, and observed that their numerical results for f(v, t) at large time could 
be collapsed on a scaling form f{v,t) ~ (v (t))~ d f(v/v (t)). Here v 2 (t) ~ (v 2 ) ~ exp[— jt], where the decay rate 
7 ~ (1 — a 2 ), represents the typical decay of the average kinetic energy. Moreover, these scaling solutions showed 
heavily overpopulated power law tails with a cross-over time, Tt a ii(f), that increases with the energy, and a is the 
coefficient of normal restitution. 

Soon after that Krapivsky and Ben-Nairn [^3| and the present authors [p4| , using a self-consistent method, gave 
a theoretical explanation of these power law tails for general dimensionality d together with explicit predictions of 
the tail exponent a, defined through /(c) ~ l/c 2a+d . These results were recently extended p5| , p6[ to driven inelastic 
Maxwell models, where the high energy tail is of exponential form exp[— A|c|]. Apparently, the type of overpopulation 
depends sensitively on the microscopic model, on the degree of inelasticity and on the possible mode of energy supply 
to the dissipative system. 

From the point of view of kinetic theory the intriguing question is, what is the generic feature causing overpopulation 
of high energy tails in systems with inelastic particles, rather than what is the specific shape of the tail. How does the 
overpopulation depend on the underlying microscopic model, and on the different forms of energy input ]^,|5|,p^-^6[ ? 
As discussed more extensively in Ref. f2q| the differences in shape are frequently related to non-uniformities in the 
limits of long times, large velocities and vanishing inelasticity, which lead to different results when taking the limits 
in different order or when taking coupled limits such as the scaling limit (e.g. the differences between bulk and tail 
behavior), or performing an expansion in powers of the inelasticity, and then studying large times (typically Gaussian 
tails are observed H]) or studying large times at fixed inelasticity and taking large time limits afterwards (typically 
overpopulated tails are observed |^,^|) with a whole wealth of coupled limits in between. 

What are inelastic Maxwell models? The IMM's, introduced in JlJJl^], share with elastic Maxwell molecules 
the property that the collision rate in the Boltzmann equation is independent of the relative kinetic energy of the 
colliding pair. However, these IMM's do not describe real particles, but only pseudo-particles with 'collision rules' 
(vi,V2) — > ( v ii v 2) between pre- and post-collision velocities, defined to be the same as for smooth inelastic hard 
spheres with a restitution coefficient a with (0 < a < 1). There are no such objects as 'inelastic Maxwell particles' 
that interact according to a given force law and that can be studied by molecular dynamics simulations. 

These IMM's are of interest for granular fluids in spatially homogeneous states, not because they can claim to be 
more realistic than IHS's, but because of the mathematical simplifications resulting from an energy-independent colli- 
sion rate. Nevertheless the IMM's keep the qualitatively correct structure and properties of the nonlinear macroscopic 
equations |27| and obey Haff's law pq] , just like the even simpler inelastic BGK or single relaxation time models p£(] 
do. What harmonic oscillators are for quantum mechanics, and dumb-bells for polymer physics, that is what elastic 
and inelastic Maxwell models are for kinetic theory. 

From the point of view of nonequilibrium (steady) states, the structure of velocity distributions in dissipative 
systems, including the high energy tail, is a subject of continuing research, as the universality of the Gibbs' state of 
thermal equilibrium is lacking outside thermal equilibrium, and a possible classification of generic structures would 
be of great interest in many fields of non-equilibrium statistical mechanics. 

After this explanation of the possible relevance of inelastic Maxwell models for different fields of research, we 
concentrate on the kinetic theory for these models, in particular on the simplest case, the freely cooling one without 



energy input. We return now to the recent results of Refs. [|21 30 22J, and we consider their observations, described 
above, as strong evidence for the existence of interesting limiting behavior when coupled limits are taken, and more 
explicitly, we interpret their findings as follows: the transformed or rescaled distribution function /(c, i) defined 
through, 

f(v,t) = (v (t))- d f(v/v (t),t), (LI) 
approaches a scaling or similarity form in the coupled limit as t — > oo and v — > 0, with v/v (t) — c kept constant, i.e. 

lim f(c, t) = hm (v (t)) d f(v (t)c, t) = /(c). (1.2) 



The coupled limit considered in (1.2) is the same scaling limit as considered in the solutions of the nonlinear Boltzmann 



equation for IHS in Refs. WJ2q]> although it has not been pointed out so emphatically. On the other hand the scaling 
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limit, considered in Refs. |13|-|l6[, is the coupled limit as t — > oo and a — > 1 with (1 — a )t = const, while u is kept 
constant. 

The present paper leads to the formulation of a conjecture by combining the recent results on scaling solutions and 
overpopulated high energy tails in inelastic hard sphere fluids and inelastic Maxwell models with an older conjecture 



of Krook and Wu |19| on the role of a special self-similar solution, the so-called BKW-mode ]17|-|20[], named after 
Bobylev, Krook and Wu. This conjecture reads: Solutions f(v,t) of the nonlinear Boltzmann equation for dissipative 
systems - i.e. the rescaled distribution f{c,t) in ( |lT| ) rescaled with the instantaneous r.m.s. velocity v (t) ~ (( l ' 2 )) 1 ^ 2 



approach for general initial conditions, in the scaling limit (1.2), to the scaling solution /(c) with an overpopulated 
high energy tail. In taking the scaling limit the degree of inelasticity, (1 — a 2 ), must be kept constant, and cannot be 
interchanged with the elastic limit (a j 1). 

This conjecture is a variation on the Krook- Wu conjecture for elastic Maxwell molecules, formulated as: "An 
arbitrary initial state tends first to relax to a state, characterized by the BKW-mode. The subsequent relaxation 
is essentially represented the BKW-mode with an appropriate phase". As it turned out, this conjecture was not 
supported by numerical and analytic results obtained for the physically most relevant initial distributions with a 
finite second moment in the limit as t — > +oo [ pof . However, Bobylev and Cercignani have recently shown that the 
conjecture holds in systems of elastic Maxwell molecules in the limit t — ► — oo for the so-called eternal solutions f(v, t) 
pl| , p2[ , which are characterized by a divergent second moment. 

The paper is organized as follows. In section II the mathematical model Boltzmann equation for IMM's is con- 
structed starting from the Enskog-Boltzmann equation for IHS's, and some basic properties are derived there, as well 
as in Appendix A. In section III we show that the IMM Boltzmann equation admits a similarity solution with a 
power law tail ~ l/c 2a+d , where the non-integer tail exponent a is the solution of a transcendental equation, which 
is solved numerically. The moments (c 2n ) of the scaling form /(c) with n < a are calculated in section IV. B from a 
recursion relation. The moments with n > a are divergent. In section IV. B we demonstrate that the moment fi n (t) 
of the rescaled distribution /(c, t), for the general class of initial conditions with all moments (v 2n ) < oo, approach in 
the long time limit for n < a to the unique set of moments ji n of the scaling form, and we analyze how all moments 
A*n(i) — ► oo for n > a. Some generalizations of these results are described in section IV. C. In section V we present our 
conclusions, and interpret our results as a demonstration of our conjecture for IMM's, i.e. as a weak form of approach 
of an arbitrary rescaled distribution /(c, i) to a universal scaling form. 

II. KINETIC EQUATIONS FOR DISSIPATIVE SYSTEMS 
A. Inelastic Hard Spheres 

For the construction of inelastic Maxwell models, it is convenient to start from the spatially homogeneous Boltzmann 
equation for inelastic hard spheres. We study the velocity distribution, /(v,r) in the so-called homogeneous cooling 
state (HCS). Here r is the "external" laboratory time, and the relation to time t used in section I will be given in due 
time. Moreover we restrict ourselves to isotropic distributions f(v,r) with v = |v| with isotropic initial conditions, 
f(v, 0). The most basic and most frequently used model for dissipative systems with short range hard core repulsion 
is the Enskog-Boltzmann equation for inelastic hard spheres in d— dimensions 0, 

d T h = 1(f) = 2 J J dv 2 6(v 12 ■ n)|v 12 • n|[^/r/r - fih], (HI) 

where /** is short for /(vj*,r), and we have absorbed constant factors in the time scale. Velocities and time have 
been dimensionalized in terms of the width and the mean free time of the initial distribution, and |vi2 • n| is essentially 
the dimensionless collision rate. Moreover, J (• • •) = (l/Clj) J dn(- • •) is an angular average over a d-dimensional unit 

sphere, restricted to the hemisphere, V12 ■ n > 0, through the unit step function 6(x), and Qd — 27r d / 2 /T(\d). 

The velocities v** with i,j = {1, 2} denote the restituting velocities, and v* the corresponding direct postcollision 
velocities. They are defined as, 

v** =v. t -i(l + i)vy -nn 

v* =v. t -i(l + a)vy nn. (II.2) 

Here a is the coefficient of restitution (0 < a < I), the relative velocity is Vy = v, — Vj, and n is a unit vector 
along the line of centers of the interacting particles. In one dimension the angular average J , as well as the dyadic 



product nn can be replaced by the number 1. One of the factors (1/cc) in Eq.(H.I) originates from the Jacobian 
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dv\*dv2* — (l/a)dvidv2, and the other one from the collision rate of the restituting colli sion s, |v£jj -n| = (l/a)|vi2 



Furthermore, in the HCS symmetrization over n and — n allows us to replace 2Q(x) in (III) by 1. 

In subsequent sections we will also need the rate equations for the average {i]j) t — J dvip(v)f(v, t), as follows from 
the Boltzmann equation, 

d(i>) T /dT = fdvi>(v)I(f) =/ n /dv 1 dv a |v ia -n|/ 1 / 2 [V(vi)-V(v 1 )]. (II.3) 

The Boltzmann collision operator conserves the number of particles (ip(~v) = 1) and momentum (V'(v) = v), but not 
the energy (tp('v) — v 2 ). Here normalizations are chosen such that, 

(1) T =/dv/(v,r) =1 
<v) r =/rfvv/(v,r) = 

<t; 2 ) T = /dvt- 2 /(v,T)=id« 2 (r). (II.4) 

As a consequence of the inelasticity an amount of energy, j(l — a 2 )[vi2 • n] 2 , is lost in every inelastic collision. 
Consequently the average kinetic energy or granular temperature (v 2 ) keeps decreasing at a rate proportional to the 
inelasticity (1 — a 2 ). So, the solution of the Boltzmann equation does not reach thermal equilibrium, described by 
the Maxwellian v?o(v) = ~K~ d l 2 exp[— v 2 ], but is approaching a Dirac delta function 8^ d \v) for large times. As the 
convergence of f[v, r) to its limiting value as r — * oo is in general non-uniform, the rescaled distribution may approach 



a different limit (see (t.2)). Also, the detailed balance condition is violated, and the Boltzmann equation does not 
obey an H— theorem. The moment equations and the behavior of the scaling solutions for freely evolving and driven 
IHS fluids have been extensively discussed both in the bulk of the thermal distribution, as well as in the high energy 
tails 



B. Inelastic Maxwell Models 



One of the difficulties in solving the nonlinear Boltzmann equation ( II. 1 ) for hard spheres is that the collision rate 



V12 • n| is not a constant, b ut is proportional to the relative velocity vi2 of the colliding pair, which is typically of order 



V12 ~ v d (t), as defined in (II.4). Maxwell models on the other hand are defined to have a collision rate independent 
of the relative energy of the colliding particles. 

In the recent literature two different types of mathematical simplifications have been introduced which convert the 
IHS- Boltzmann equation into one for an inelastic Maxwell model with an energy independent collision rate. In the 
most drastic simplification, the IMM-A discussed in Rcfs. |3(i| , p3| , p4| , one replaces the collision rate |via • n| for the 
direct collisions, as well as the one for the restituting collisions, |v*2 ■ n| = |vx2 1 n |/ a by ^ s typical mean value v (t). 
In a more refined approximation Bobylev et al. [jl3| replace these collision rates by u (r)|vi2 • n|. We call this model 
IMM-B. Both approximations keep the qualitatively correct dependence of the total energy v 2 (t) on the 'external' 
time r. In fact, strictly speaking these models should be called pseudo-Maxwell molecules, because there do not exist 
microscopic particles with dissipative interparticle forces, for which the mathematical model Boltzmann equations 
below can be derived. 



By making the above mathematical simplifications we obtain from (III) a collision term which is multiplied by a 
factor v {t). This factor is then absorbed by introducing a new time variable t. For model IMM-A the resulting time 
transformation and Boltzmann equation in dimensionless variables are then given by, 



dt = v {r)dT 



dth = i(f) 

= -h + J J rfv 2 i/r/r- (n.5) 

For the IMM-B we introduce a slightly different time variable t, and obtain the Boltzmann equation in dimensionless 
variables, 

dt = f3iv (r)dT 

dth = i(f) = J J rfv 2 | Vl2 • n | [i/r/r - hh] 

= -/i+/ /dv 2 |v 12 -n|i/r/r, (H.6) 
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where f' = (l//?i) / with /3i defined in (A.l). The prefactors in the time transformations are chosen such that the 
loss term takes the simple form —ft. This implies that the new time variable t(r) counts the average number of 
collisions suffered by a particle within the "external" time r. Hence t is the collision counter or "internal" time of a 
particle. 

One of the important properties of Maxwell models is that the moment equations form a set of coupled equations, 
that can be solved sequentially. For the one-dimensional Maxwell model these equations were derived in Ref. [fl2|| , 
and for the three-dimensional model IMM-B with uniform impact parameter this was done in Ref. |l3| . The general 
moment equations for the present Maxwell models will be derived after having obtained the characteristic function 
in Sect. III. At this point we make an exception for the second moment, which determines the typical velocity v a (t) 
through the relation (v 2 ) t = ^dv 2 (t), needed to study the rescaled distribution function. It follows from the kinetic 
equations (II. 5) and Qll.6|) as 



d t {v 2 )r- 



~l(v 2 )t 

' l-a 2 . 
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2p(l~P) 
d 



2d 

l-a 2 _ 4p(l-p) 



d+t 



d+1 



(IMM-A) 
(IMM-B) 



(II.7) 



where p = \(1 + a). Consequently v a (t) = v o (0) exp[— hjt]. Moreover, by solving the differential equations for t in 
(|II.5|) and ([1. 6) we obtain the relations between the internal time t and the external time t, i.e. 



exp[±7t] 



±7v o (0)t (IMM-A) 
i/3i 7 t> o (0)r (IMM-B) ' 

as well as the decay of the energy in terms of internal time t and external time r, i.e. 




v 2 {t) = exp[-7t]i£(0) = 



v 2 o(0)/[l + h v o(V)T] 2 (IMM-A) 
^(0)/[1 + ±/3 i7 i; o (0)t] 2 (IMM-B). 



(II. 



(II.9) 



This shows that Haff's law pq] , given by the second equality, is also valid for inelastic Maxwell models. 

To further elucidate the difference between the two classes of models we change the integration variables n — which 
specifies the point of incidence on a d-dimensional action sphere of two colliding particles — to the impact parameter, 
b = |vi2 x n| = sin6>, where 9 = cos -1 (vi2 • n) is the angle of incidence. The relevant structure of the integrals for 
IMM-A and IMM-B is respectively, 



V12 • n 



d0(sin0) d - 2 ~ / dbb d - 2 /^l-b 2 
(t Jo 

tt/2 f l 

d9 (sin 9) d ~ 2 cos 9 ~ / dbb d ~ 2 



(11.10) 



Therefore model A has a uniform distribution V(9) — 1 over angles of incidence, and a non-uniform distribution 
P(b) = 1 / V 1 — b 2 over impact parameters, which is biased towards grazing collisions, where 6 = 1. Model B has a 
uniform distribution P(b) = 1 over impact parameters, and a non-uniform distribution, V(0) = cos 9, biased towards 
zero angle of incidence. In this context it is important to note that the arguments on the validity of the Boltzmann 
equation are based on the assumption of molecular chaos, i.e. absence of precollision correlations between the velocities 
of a colliding pair. This implies that the distribution function of impact parameters be uniform. Hence, model IMM-A 
does not obey molecular chaos. 

The question of interest in then: do IMM-A and IMM-B yield qualitatively the same results for the scaling distribu- 
tion? The question is rele vant because Molecular Dynamics simulations (MD) of a system of N IHS have shown that 
the dissipative dynamics (II. 2) drives an initially uniform distribution P(b) = 1 in the HCS towards a non-uniform 
distribution P(b), biased towards grazing collisions, which violates molecular chaos. So the IMM-A model with a 
built-in initial bias may lead to spurious effects, such as power law tails in f(v, t), which are artifacts of a too drastic 
simplification. 



C. Similarity Solutions 

The questions addressed in this subsection are: do the Boltzmann equations for the Maxwell models, constructed in 
the previous section, admit similarity solutions, and what are the properties of such solutions? We define a similarity 
solution /(c) through the relation, 
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f(y,t)=v- d (t)f(v/v (t)). 



(11.11) 



The normalizations imposed by (II. 4 ) on these solutions are, 



dc/(c) 



dc c 2 /(c) 



¥■ 



(11.12) 



By inserting (11.11) in (II. 1), and using v Q (t) ~ exp[— hjt] we obtain the following integral equation for /(c), i.e. 



I 7 1 

2 '9c 



c/(c) = /(/). 



(11.13) 



Here the operator /(/) has the same functional form as /(/) in ( JII.5| ) or (II.6) with {vj, /} replaced by {c», /}. 

One of the goals of this paper is also to analyze in section IV in what sense the rescaled distribution function, 
f(c,t), approaches its limiting form as t — > oo. To do so, we also need the kinetic equation for the rescaled f(c,t), 
which reads, 



ftZ+M-c/ = /(/)• 



(11.14) 



Some comments are in order here. Physical solutions /(c) of (11.13) must be non-negative. A velocity distribution 



/(y, t) , evolving under the nonlinear Boltzmann equation, pr eserve s positivity for a positive initial distribution /(v, 0) 
|53],[34j]. However, for scaling solutions, being the solution of ( 11.13 ), positivity is not guaranteed |]l3| , p0| . If one would 
know a positive scaling solution /(c) - as is the case in on e dimension |H| - and prepare the system in this initial 



state, then the entropy S or the H— function in state (11.11), also shows singular behavior, i.e 



S(t) = -H(t) = - J rfc/(c) In /(c) - \d^t + const 
~ —7}djt + const (t large), 



(11.15) 



where 7 is positive and J /In / is some constant. In these solutions the entropy keeps decreasing at a constant rate 
id 7. This is typical for pattern forming mechanisms in configuration space, where spatial order or correlations are 
building up, as well as in dynamical systems and chaos theory, where the rate of irreversible entropy production is 



negative on an attractor pq-|37|. In fact the forward dissipative dynamics, (v l5 v 2 ) — > ( v i\ v 2)i defined in (II. 2), has 
a Jacobian J = a < 1, i.e. dv^dv^ — adv\dv2 1 corresponding to a contracting flow in v— space. Moreover, there 
is no fundamental objection against decreasing entropies in an open subsystem, here the inelastic Maxwell particles, 
interacting with a reservoir. The reservoir is here the sink, formed by the dissipative collisions, causing the probability 
to contract onto an attractor. 



III. POWER LAW TAILS 



A. Fourier transformed Boltzmann equation 

The goal of this section is to show that the Boltzmann equation for IMM's has a scaling solution with a power law 
tail. This is done by introducing the Fourier transform of the distribution function, (p(\s.,t) = (exp[— ik • v]) t , which 
is the characteristic function or generating function of the velocity moments. Because f(v,t) is isotropic, ip(k,t) is 
isotropic as well. It is also convenient to consider <j){x, t), defined through the relation (p(k, t) = <j>(jk 2 , t). 

W e start with the simplest case, and apply Bobylev's Fourier transform method fll3[fi1j| ] to the Boltzmann equation 
(fL5|) for model IMM-A with the result, 



ftp(k,t)= [ Kk+,<Mk_,t)-^(0,iMk,i)] 

J n 

dt<t>(x,t)= I [<l>(xe + (n),t)<t>(xe-(n),t) - ct>{0,t)<l>(x,t)}, 



(III.1) 



where <p(0,t) = 1 = cf>(0,t). Here we have used ( II. 3 ) with ^(vi) = exp[— ik • v*] and expressed the exponent as 
k • v* = k_ • vi + k + • v 2 (see (tl.2)), where 



G 



fck_, 



pk ■ nn 
k k. 



k\ =p 2 (k-n) 2 =e+(n) 
fc 2 = [1 - z(k- n) 2 ] = e_(n), 



k = fck_ 

with p — ^(1 + a) and z — 2p — p 2 . In one dimension this equation simplifies to 

d t (p(k,t) = ip{pk,t)tp({l-p)k,t) - <p{k,t), 



(HL2) 



(III.3) 



where fc + = pk and fc_ = (1 — p)fc. Equation (ill.l) has the interesting property that for a given solution ip{k,t) 
one has a whole class of solutions </>(k, t) = exp[ik • w]<p(k, i) where w is an arbitrary velocity vector |18 13}| , This 
property reflects the Galilean invariance of the Boltzmann equation. 

Because f(y, t) is isotropic, only its even moments are non- vanishing, and the moment expansion of the characteristic 
function then takes the form, 



v(M)= E 



(-ife) T 



m— even 



nil 



.((k. v r) t = E 



-m„(i) = (j){x,t), 



(III.4) 



where x = fc 2 /4. The angular time independent average ((k-v) 
is defined as 



2r>\ 



/?2n is calculated in (A.l) and the moment m n (t) 



m n (t) = 4 n nl(3 2n (v 2n ) t /(2n)l = (v 2n ) t /(d/2) n 
(3 2n = ((k ■ v) 2 "> = (l/2) n /(d/2)„ 



(III.5) 



The Pochhammer symbol (a)„ is defined in (A. 3) and we have used the duplication formula for the Gamma function 
(2n)! = r(2n+ 1). Furthermore we note that the moments of a Gaussian ip (v) = ir~ d l 2 exp(— v 2 ) are (v 2n ) Q — (d/2) n . 

Scaling solutions in Fourier representation have the form (f)(x, t) = $>(e (t)x), where $(-|fc 2 ) is the Fourier transform 
of /(c) and e a (t) = v 2 (t) — e rt e o (0). Substitution in Eq.( III.l ) yields the integral equation for the scaling form, 



■ jx&(x) + <f>(x) = / $(xe + )<i>(xe_) 

J n 



(HI.6) 



Its moments follow from the expansion, 

[in = Jdcc 2n f(c)/(d/2) n = (c 2n )/(d/2) r 

where ( [1.12) imposes fix = 1. 



(HL7) 



B. Small-fc singularity of characteristic function 



In case all coefficients in the Taylor expansion (III. 7) exist, then $(x) is regular at the origin, and the corresponding 
scaling form /(c) falls off exponentially fast at large c, and all its moments (c 2n ) are finite. Suppose now that the 
small-fc or small- a: behavior of $(x) contains a singular term x a , where a does not take integer values (note that even 
powers fc 2a = fc 2n represent contributions that are regular at small fc), then its inverse Fourier transform scales as 
l/c 2a+d at large c. For this distribution the moments with n > a are divergent, and so is the n-th derivative of the 
generating function <E>(x) at x — 0. The requir emen t that the total energy be finite imposes the lower bound a > 1 



on the exponent because of the normalization (11.12) 



Model IMM-A 



To test whether the Fourier transformed Boltzmann equation (III.l) admits a scaling solution with a dominant 
small-x singularity x a , we make for the small— cc ansatz, 



1 



Ax a , 



(J1U 



insert this expression in ([II. 6), and investigate whether the resulting equation admits a solution for the exponent a. 
This is done by equating the coefficients of equal powers of x s on both sides of the equation, which yields for general 
dimensionality, 



7 



7 = Ai= /"[l-e+(n)-e_(n)]=2p(l-p)/d 

J n 



07 : 



[l- e+ (n)-e«(n)]. 



(III.9) 



The eigenvalue A s (s = 1, o) has been calculated in ( A. 4 ) and flA.q ) of the Appendix. The relevant properties are: (i) 
Ao = because of particle conservation; (ii) lim s ^o A s = — 1; (Hi) A s is a concave function, monotonically increasing 
with s, and (iv) all eigenvalues for non-negative integers n are positive (see Fig.l). In one dimension the above 
eigenvalue becomes A a = 1 — p 2a — (1 — p) 2a . 




FIG. 1. Eigenvalue A s is a concave function of s, plotted for different values of a for the 2-D inelastic Maxwell model IMM-A. 
The line y = sAi is plotted for a = 0.6, 0.8 and a — l(top to bottom). The intersections with A s determine the points so (filled 
circles) and si (open circles). Here Si = o determines the exponent of the power law tail. For the elastic case (a = 1, Ai = 0, 
energy conservation) there is only I intersection point. The branch of A s in the interval so < s < si is referred to as stable and 
the branches s < so and s > si as unstable. 



To continue we combine both relations in (III. 9), which determine the exponent a as the root of the transcendental 
equation, 

A s = sAi. (III. 10) 

This equation has been solved numerically, and the results are plotted in Fig. 2. We note here that Krapivsky and 
Ben-Nairn p3| have derived the same transcendental equation. 



As can be seen from the graphical solution in Fig.l, the transcendental equation (111.10) has two solutions, the 
trivial one (s = 1) and the solution s\ = a with a > 1. The numerical solutions for d = 2, 3 are shown in Fig. 2a as 
a function of a, and the a-dependence of the root o(a) can be understood from the graphical solution in Fig.l. In 
the elastic limit as a | 1 the eigenvalue Ai(a) — > because of energy conservation. In that limit the transcendental 
equation ( IITlOj) , A s (l) — sAi(l) = 0, no longer has a solution with a > 1, and a(a) — ► oo, as it should be. This is 



consistent with a Maxwellian tail distribution in the clastic case. Needless to say that the transcendental equation 
for the one-dimensional IMM-A has the solutions sq = 1 (trivial) and s\ = a = 3/2, describing a power law tail 
/(c) ~ _L/\c\ 2a+d = |c| 4 , in full agreement with the exact scaling solution /(c) = (2/7r)(l + c 2 ) -2 , found by Baldassarri 
ct al. [El] for this case. The one-dimensional case is a bit pathological because the intersection points (sq, si) = (1, 3/2) 
of y = A s (a) and y = sAi(a) are independent of p = 5(1 + 0;), implying that these points are common to all A— curves 
at different parameter values a. In summary, we conclude that there exists for inelastic Maxwell models a scaling 
solution /(c) with a power law tail l/|c| 2o+d at large energies. 

As a parenthesis we apply the previous analysis for similarity solutions to the elastic case (a = 1), where the BKW 



- mode, discussed in the introduction, is an exact similarity solution. In the elastic case similarity solutions of ( III. 1 ) 
would also have the form <p(x,t) = ^(e (t)x), but the time dependent factor e (t) is not determined by the conserved 
second moment (v 2 ) = |d So the typical time scale 7, entering through e {t) = exp(— jt), will be determined by 
the lowest moment with t— dependence, i.e. (i; 4 ) ~ £ 2 (t), and its rate equation imposes A2 = 2j. The transcendental 
equation becomes then 7^ = A s — SA2/2, which has again two solutions, sq = 2 and si = 3 because A2/2 = A3/3. 
As both a-values are integers, the small— a; behavior of the characteristic function contains only regular terms x 2 and 
x 3 , which do not result in any power law tails. In fact, the solutions {s = 2,3} correspond to the exact closed form 
solution (f>{x) = e~ x (l + x), the well-known Bobylev-Krook-Wu mode for elastic Maxwell molecules Jl8|,p0|. 
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FIG. 2. Exponent a(a), which is the root of ( HI.1C ), as a function of the coefficient of restitution a, which determines the 
high energy tail l/c 2a+d of the scaling solution /(c). Left panel for IMM-A with uniform angle of incidence in 2-and 3-D, and 
right panel for IMM-B with uniform impact parameter in 3-D. 



Model IMM-B 



The second part of this section deals with the Boltzmann equation (H.6) for model IMM-B with uniform impact 
parameters, as introduced by Bobylev et al. |l3fl , and we show that the above method gives similarity solutions with 
power law tails for this model as well. We restrict the analysis to the three-dimensional case. The integral equation 
for the characteristic function in this case has been derived in Et], and reads, 



6W(k, t) 



|k.n|0(k + ,t)0(k_,i)-0(k,i). 



(HI.11) 



In a similar way as in the previous section, we derive the integral equation for the corresponding scaling solution, 
4>{x,t) = $>(xe (t)) with the result, 



- 7X$'(x) + $(x) = / |k-n|$(xe + )$(ze_), (111.12) 

J n 

and the procedure described in ([II. 8) and (III. 6) leads to the same transcendental equation ( IILlCj ) with eigenvalue, 



Ik-nlfl-e! 



2 ( dxx{l-p 2a x 2a - [1 - zx 2 ] a } 
Jo 

i ri- (i-p) 2a + 2 



= i - 



a + 1 a + 1 [ 1 - (1 - p) 



(III. 13) 



The equality on the second line has been obtained by changing to the new integration variable x — cos 9, and the 
value for A a is in agreement with the c orres ponding eigenvalue obtained in |]l3|| for integer a. For a = 1 one obtain s 
here Ai = p(l — p), in agreement with (II. 7) for d=3. The numerical solution of the transcendental equation ( III. 1C ) 
is shown in Fig. 2b. 

It turns out that the numerical values a(a) for the exponents of IMM-B with a uniform distribution of impact 
parameters looks qualitatively the same as those for model IMM-A with a uniform distribution of angles of incidence. 
The differences in distribution of impact parameters of both models, uniform in IMM-B versus biased towards grazing 
in IMM-A, has no qualitative effect on the nature of the singularity, i.e. on the values or a-dependence of the tail 
exponents. Therefore, the power law tail in /(c) ~ l/c 2a+d is not a spurious effect induced by models with impact 
parameters b biased towards grazing collisions. 

Afte r completion of this article Bobylev and Cercignani p8| have kindly informed us that equations (III. 11) and 
(111.12) are also correct for the d— dimensional version of model IMM-B. This implies not only that the exponent a 



<) 



can be calculated for general dimensionality, but it also implies that the next section, after some trivial modifications, 
applies to the model IMM-B in d dimensions as well. 

We conclude this section by discussing some numerical evidence for the conjecture on the approach to a scaling 
solution with algebraic high energy tails, as formulated in section I. Baldassarri [0 has obtained large-i solutions 
f(v,t) by applying the DSMC (Direct Simulation Monte Carlo) method to the Boltzmann equations for three types 
of inelastic Maxwell models, among which the two-dimensional IMM-A model and the three-dimensional IMM-B 
model, analyzed in this article. For the totally inelastic case (a = 0) he has observed that f(v,t), for general initial 
data, evolves after sufficiently long time to a scaling solution of the form ( |l.l[ ), on which the simulation data can be 
collapsed. Moreover it has a power law tail, /(c) ~ l/c 2a+d in agreement with the predictions in Figs. 2a, b at a = 0. 



IV. APPROACH TO SCALING SOLUTIONS 



A. Moment Equations 



The moment equations for Maxwell models are special because they form a closed set of equations that can be 
solved sequentially as an initial value problem. In this section we study the effects of power law tails /(c) ~ A/c 2a+d 
on the moments, and investigate in wh at sense, if any, the calculated time dependence of the moments as t — > oo is 
related to the singular behavior ([II. 8) of the scaling form $(x), derived in section III. The latter implies that the 
moments [x n , generated by $(x), with n > a are divergent and that those with n < a remain finite. 

First consider the standard moments m n (t). Inserting the expansion (III. 4) in ( [II. 1 ) and equating the coefficients 
of equal powers of x yields for the moments m n (t) 



,,2n\ 



■ /{d/2) n the following equations of motion, 



777, 



+ A„?77„ = ^-ff (Z, 77 - /)777;777 r . 



(IV.l) 



/=1 



where the coefficients H(l,m) and eigenvalues A„ for model IMM-A are defined and calculated in ( A.4)-( [A.8[ ). Those 
for model IMM-B fo llow after some trivial replacement, indicated in Appendix A. Re gard ing the moments, we have 
777 = 1 because of ( pi. 4| ), and m\{t) — exp(— Ait)m 1 (0), where Ai = 7, as given in (II. 7). Moreover as t — > 00, all 
moments with 77 > vanish, which is consistent with the limiting behavior f(v,t) — > 5 W (v) for t — ► 00. 

Next we consider the moments fi n generated by the scali ng form <j>(x). Using a self-consistent argument we have 
demonstrated in section III that the kinetic equation (II. 5) admits a scaling solution $(x) with a dominant small- 
x singularity x a , with a > 1 and non-integer. This implies that all 77— th order derivatives of $(x) at x = 0, or 
equivalently all moments fj, n , are finite if n < n a = [a] < a, and all those with n > a are divergent. Here [a] is the 
largest integer less than a. Hence, the small-x behavior of $(x) can be represented as , 



$(x) = £(- 

n=0 



\n Mr. 



o(x a ), 



(IV.2) 



where the remainder is of order o(x a ) as x — > 0. In this scaling form we only know the exponent a and the moments 
Mo = Ml = !• Now we calculate the unknown finite moments of the scaling form, /j„ with 1 < n < a. This is done by 
inserting (IV.2) into the kinetic equation ( III. 1 ) , yielding the recursion relation, 



1 = 1 



7n = A„ - 777 = A„ - 77Ai. 



(IV.3) 



Here 71 = 7 — Ai = on account of (II. 7) and the initialization is /Ji = 1. The solutions pi n for 77 = 2, 3, 4, 5 in model 
IMM-A are shown in Fig. 3 as a functio n of the coefficient of restitution a. Furthermore we observe that the root s = a 
of the transcendental equation ( III.lQj ), 7$ — A s — sAi = 0, indicates that 7^ changes sign at si = a (see open circles 
in Fig.l), and that according to Section III all moments fx n with n > s\ — a (n on unstable branch) are divergent. 
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FIG. 3. The moments fi n ~ (c 2n ) of the scaling form /(c) for n = 2, 3, 4, 5 as a function of a. The moment /X2 exists for all 
values of a, while higher moments only exist for a above a threshold value, indicated by the vertical asymptotes. 



The recursion relation ( tV.3 ) for the moments fi n in the one-dimensional case is again a bit pathological in the sense 
that the stable branch (s = 1 < s < Si = 3/2) contains only one single integer label, i.e. s — 1. So only fi = /i X = 1 
are finite, and all other momen ts are infinite, in agreement with the exact solution of Baldassarri et al. 

The recursion relation (1V.3) for both models has a second set of solutions which has been studied by Bobylev 

ct al. filfl, w ho showed that within this set there are moments fi n being negative. The argument is simple. Consider 
fi n in (IV. 3) with n = n + 1. Then the prefactor l/j no +i on the right hand side of this equation is negative because 
the label n a + 1 > a is on the unstable branch of the eigenvalue spectrum in Fig.l, while all other factors are positive. 
This implies that the corresponding scaling form /(c) has negative parts, and is therefore physically not acceptable. 
We also note that the moments fi n of the physical and the unphysical scaling solution $(x) coincide as long as both 
are finite and positive in the a— interval that includes a = 1. 



B. Long time behavior of rescaled moments 

In the introduction we have already indicated that the distribution function f(v,t) itself approaches a Dirac delta 
function, <5^ d ^(v), and that the rescaled distribution function f(c,t) — v^(t)f(v (t)c,t), as denned in (|L"T|), supposedly 
approaches a scaling form. To demonstrate how this happens, we define the rescaled function 4>(x, t) through the 
relation 4>(x,t) — <p(e (t)x,t) with e (t) = exp[— jt]e o (0). We also consider the more general case where the typical 



time scale is not a priori fixed by defining e a (t) in terms of one specific moment (v ). Then it follows from ( III. 1 
that 

d t 4>(x,t)-^xd x 4>(x,t)+4>(x,t)= 4>(xe + ,t)<t>(xe-,t), (IV.4) 

J n 

where 0(x,t) with x — jk 2 is the Fourier transform of f(c,t) in ([Ll]). The stationary solution of this equation 
determines a possible similarity solution for a given 7. The equations of motion for the rescaled moments jJL n (t) of 



f(c,t) are obtained by substituting the Taylor expansion 4>(x,t) = *Yl, n (—x) n iJL n (t)/n\, into ( IV.4 ), which yields for 

71=1,2,... 



n-l 



fin + Jntln = ^ H ( l > 71 ~ OWMn-J 



2=1 



7„ = A„ - 717, (IV.5) 

with /j„(0) = m n (0)/e™(0). We also observe that the equations for <j>(x, t) and for the moments fi n (t) hold for arbitrary 
rescaled initial distributions, /(c, 0) = v% (0)f(cv o (0), 0). 
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FIG. 4. Graphical construction of the solutions of A a = S7 for different values of 7 at dissipation a — 0.5. Intersections (so) 
with smaller s— values are marked with filled circles, and intersections (si) with larger s— values with open ones. 



To study the long time approach of f(c,t) to the scaling solution /(c) in (|Ll"|), we analyze the approach of the 
rescaled moments fi n (t) = m n (t)/e"(t) to the moments fi n of the scaling solution f(c). To do so we have to choose 



7 = Ai (see Fig. 4) on account of (II. 4) and (II. 7). The infinite set of moment equations (IV. 5) for /j, n (t) can be solved 
sequentially for all n as an initial value problem. To explain what is happening, it is instructive to use a graphical 
method to determine the zeros of the eigenvalue 7^ = \ s — sj for different values of 7. This is illustrated in Fig. 4 by 
determining the intersections {sq, si} of the curve y = X s and the line y = sj, where sq and s± are denoted respectively 
by filled (•) and open circles (o). These circles divide the spectrum into a (linearly) stable branch (so < s < s±) and 
two unstable branches (s < sq) and (s > s±). The moments /J, s (t) with s = n > a are on an unstable branch (<y s < 0) 
and will gr ow for large t at an exponential rate, fi n (t) ~ /x n (0) exp[|7„|t], as can be shown by complete induction 
from ( jTVlg ) starting at n = [a] + 1. They remain positive and finite for finite time t, but approach +00 as t — ► 00, in 
agreement with the predictions of the self consistent method of Sect. III. The moments with n on the stable branch 
are linearly stable (7^ > 0), but may still grow through nonlinear couplings with lower moments ni(t) whenever 
I = 1, 2, ■ ■ ■ is on the unstable branch (I < sq). This is relevant for the discussion in the next subsection. In the case 
7 = Ai under consideration however, all moments fii with I = 1, 2, • • • , [a] are globally stabl e and approach for t — > 00 
the limiting value /x n (oo) = fj, n , which are the finite positive moments of the scaling form (IV. 2), plotted in Fig. 3. 

The behavior of the moments described above is considered as a weak form of convergence or approach of f(c,t) 
to f(c) for t — ► 00. This result will be described in a more precise form and summarized in the last section. The 
physically most relevant distribution functions f(v, t) and rescaled distribution functions f(c, t) are those with regular 
initial conditions, i.e. all moments m n (0) = e"(0)/i„(0) < 00. This implies that the initial condition for the rescaled 
Fourier transform <p(x,0) with x = jk 2 of /(c, 0) has a series expansion regular at the origin x — 0, i.e. all its 
derivatives exist in that point. 



C. More similarity solutions 



The results of the previous subsection can also be generalized to different values of 7 and different classes of initial 
distributions, which possess only a finite set of bounded moments m„(0). Here we summarize only the most important 
results of these generalizations, the details of which will be published elsewhere p9| . 

The case, 7 > Ai, concerns the line y — sj2 in Fig. 4. Then fii(t) will always be in the unstable region (there 
is either an intersection point sq > 1, or there are no intersection points at all). So, (J>i(t) — * 00 as t — * 00 at an 
exponential rate, and it drives all moments fx n (t) — > 00 through nonlinear couplings to H\(t), which is present on the 
right hand side of (IV. 5). Consequently, there is no approach to the scaling form <J>(cc) with a rate constant 7 > Ai. 

The case, < 7 < Ai, concerns the line y — S74 where there are two intersection points {sq, s\} with so — a < 1 and 
s% = b > 1, corresponding to the dominant singularity x a and a subleading singularity x . The dominant singularity 
corresponds to a power law tail /(c) ~ l/c 2a+d with a < 1. Because a < 1, the system has infinite energy at all times. 
All higher moments \i n with n = 2,3,... are divergent as well at all times. Consequently, the initial states under 
discussion are already singular with cf>(x, 0) ~ 1 — Ax a . Of course such states are of much less interest for possible 
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physical applications than the regular ones, discussed in the previous subsection. The feature of interest here is to 
demonstrate that inelastic Boltzmann equations generate for initial states characterized by a singularity x a with a < 1 
a new type of singularity x b , which is found through the graphical construction using line y = S74. For elastic Maxwell 
molecules such states have been analyzed recently by Bobylev and Cercignani |[lj . In one-dimensional systems these 
initial states are closely related to Levy distribution £□] with the characteristic function tp(k) = cxp[— bk 2a ], where 
b is positive. For such distributions it is well known that Fourier inversion yields for < a < 1 a non-negative 
distribution f(v) with a power law tail \/ v 2a + d with d — 1. On the other hand, for a > 1 Fourier inversion may lead 
to a distribution with negative f(v). 

However, in this case one can say more. Following Bobylev and Cercignani, we assume that we can construct a non 
negative initial distribution with a cf)(x, 0) which is a regular function of £ = x a in a neighborhood of £ = 0, i.e. 

#r,0) = f>ir^L n (0), (IV.6) 

where |^4„(0)| < A", and we have chosen normalization such that A o (0) = Ai( 0) — 1. We have slightly modified the 
example of |31| to have a finite non- vanishing radius of convergence of fvlf . Then we can demonstrate that the 



rescaled characteristic function <j){y,t) = (f)(x,t) approaches in the scaling limit (t — > 00 with y = xexp[— yt] fixed) a 
positive scaling form or similarity solution <fr(y) with 7 = A a /a, and 

<%) = £(-l)»^ n (oo) + o(y»), (IV.7) 

71=0 

where Aq(oo) — A±(po) = 1. The approach is again in the weak sense of a finite set of moments. The ^4 n (oo) are 



positive and can be calculated from a set of recursion relations, rather similar to (CV.S). Moreover, b = s\ is the left 
most intersection point of X s with the line y = S75 (see the left most open circle on X s ). So, the initial 4>{x, 0), which 
is regular in £ = x a around ( = 0, develops a new singularity of the type x b . This can again be demonstrated by 



considering the rescaled function <j>(y,t), defined as <j){x,t) — <fi(e 7t y 7 t), and expanding (j>(y 7 t) in a series like (IV.6) 



with A n (0) replaced by A n (t). The coefficients satisfy moment equations, rather similar to (IV. 5). 

In the case, 7 < 0, the results are similar to those in the previous paragraph, except that there is only one 
intersection point at so = a (see the line y — 575 = — |7s|s). The energy and all higher moments are infinite, and the 
scaling form of the characteristic function <fr(x) is a regular function of ( = x a near the origin. A similar solution for 
the elastic case has been obtained in Ref. Bill. 



V. CONCLUSIONS 



Using self consistency arguments we have shown that the nonlinear Boltzmann equation for the inelastic Maxwell 
models IMM-A and IMM-B admits a scaling solution /(c) with a power law tail l/c 2a+d where the exponent a is 
given by the root a with a > 1 of a transcendental equation. This implies that all moments fx n of f(c) are divergent 
for n > a, and those fj, n with n < a are positive and finite, and are given recursively through (IV. 3) with initialization 
Mi = 1- 

For systems with dissipative dynamics we have formulated a conjecture about the long time approach, for general 
classes of initial distributions, of the rescaled distribution function /(c, t) — v^(t)f(cv (t),t) to a scaling solution f(c) 
with an high energy tail, which is overpopulated when compared to a Gaussian distribution. For the inelastic Maxwell 
models, studied in detail in this article, we have demonstrated this approach to the scaling solution /(c) in subsection 
IV. B in the following weak sense: 

Given on the one hand the rescaled moments fi n (t) = m n (t)/e"(t) of the physically most relevant class of regular 
initial distributions with all moments m n (0) bounded, and rescaled according to (|I.l|), and given on the other hand 
the unique set of finite positive moments /i„ of the scaling form with n < a, then all ii n (t) with n < a approach for 
t — ► 00 the limiting value /x ra through a sequences of positive numbers, and all moments fi n (t) with n > a behave as 
ix n {t) ~ /i n (0) exp[|7„|<] — > +00 as t — > 00. This is in agreement with all known properties of the scaling solution. 

We consider the long time behavior of the set of moments /i n (£) — * jj, n (n < a) and fi n (t) — > +00 (n > a) as a 
demonstration that f(c,t) approaches the scaling form /(c) for t — » 00. We further note that these results for the 
long time behavior of inelastic Maxwell models, after the rescaling of the initial distribution, are universal, i.e. they 
are independent of all details of the initial distributions. We consider the present results for inelastic Maxwell models 
as strong support for our more general conjecture, which is also confirmed for inelastic Maxwell models through the 
Monte Carlo simulations of the Boltzmann equation by Baldassarri et al, in which the approach of /(c, t) to a positive 
power law tail with the predicted exponent a was confirmed within reasonably small error bars. 
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APPENDIX A: ANGULAR AVERAGES IN D-DIMENSIONS 



The angular average ( III. 5 ) of powers of a • n = cos 9 can be simply calculated by using polar coordinates with a as 
the polar axis. Then, 



1% 



f£ de (sin ef- 2 (cos ey _ r(^±i)r(|) 



r(^r(|) 



(A.l) 

(A.2) 
(A.3) 



In fact, we will use the notations (a) s , /32s, s\ = (l) s and I ^ ) also for non-integer values of s by expressing these 
quantities in terms of Gamma functions. 



1 ' J*d6(sm6) d - 2 
For s = In this formula can be expressed as 

[3 2n = (l/2)„/(d/2) B , 

where the Pochhammer symbol is defined as 

(a) n — T(a + n)/T(a) — a(a + 1) . . . a(a + n — 1) 



Next we define H(l,m) and A„, introduced in (IV. 1) for model IMM-A, 

H(l,n-l)= ^ Jj l + (n)e n _- l (n) 
X„ = l + S n0 - H(0, n) - H(n, 0) = / [1 + 5 n0 - ei(n) - eH(n) 



(A.4) 



where e±(n) has been defined in ( III. 2 ). For model IMM-B one needs to replace J with J \a- n|. These expressions 
hold for d = 1,2, Evaluation of H(l, to) requires, 



/(to, n) 



O n 2m fir/2 

J _ 1 / ^(sin#) d - 2 (cos0) 2m [l-zcos 2 0]' 



(A.5) 



where z = 1 — (1 — p) 2 . Following Krapivsky and Ben-Nairn J2j| we change to the new integration variable /x = cos 2 9, 
to find for d = 1, 2, . . . 



H(m, n) — 



f n 


+ m 




m 


/ n 


+ TO 




m 


= p 2m 


c 



p 



2rn 



B(-^, i) Jo 



dw m - 1/2 {l - m)^[i - zn\ n 
1 d 



1=0 



(A.6) 



On the second and third line we have used the fundamental integral representation for the hyper-geometric function 
2-Fi, and its Gauss series, i.e. 



2 F 1 (a,b,c;z) 



B^i^c-b) dtt b - 1 {\ - f) c - fc - 1 (l - zt)- a 



(A.7) 



When a = — n, (n = 0, 1, 2 . . .), then 2-Fi( — n, 6, c; z) is a polynomial of degree n in z, and the Gauss series ends at 
/ = n. 

To calculate the eigenvalue X n in (A.4) we deduce from (A.7) 

H(n, 0) = I3 2n p 2n 



H(0,n)= 2 Fi(-n,i |;*) = £( ")^ (_2) 



;=o 



(A. 
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